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I. INTRODUCTION 



A first-principles, lattice QCD study of B° — B° mixing puts the Standard Model to a test 
and constrains the physics beyond it. In particular, it provides a number of theoretical in- 
put parameters necessary for constraining the elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. In the Standard Model, the oscillation frequency of B° — B° mixing is 
described by the following expression p|: 



Am„ 



Gprn^So 
lQn 2 m B 



\V*V tb \ 2 r] B M 



(1) 



where q G {d, s}, V tq and V t b are CKM matrix elements. Except for the matrix elements and 
M q all quantities in Eq. are known from experiments or can be accurately computed by 
analytic means. The hadronic matrix element 



M q = (B° q \ [bY(i - 75MM1 - 7 5 )g] \B° q ) 



(2) 



due to its nonperturbative nature is best determined from a lattice calculation. However, 
the systematic errors inherent to lattice calculations as well as the uncertainties associated 
with the perturbative matching between continuum and lattice operators limit the precision 
with which M q can be determined. A more promising approach would be to focus not on 
M q itself but on the following ratio 



Am s 
Am d 



m Bd Mi 
m Ba M t 



(3) 



This ratio relates the matrix element \Vtd\ to the experimentally measurable quantity 
Am s / Arrid- For historical reasons the hadronic matrix element is often factorized as 
M q = (8/3)m B f\ B Bq - Using this convention we can rewrite Eq. Q as 



Am s 
Am d 



V td 



m Bs ^ 2 



(4) 



where 




fh d V B Bd ' 



(5) 



whose value is directly used in CKM matrix analysis |2j. Note that since £ is the ratio of 
the same quantities with only the mass of the light quark being changed, many systematic 
errors in the lattice calculation should cancel. Therefore, the focus of this paper will be the 
determination of the ratios f B J f Bd , B Bs /B Bd and £, and not the individual quantities f B 
and B B . Moreover, we will calculate £ in two ways. An "indirect" approach is to calculate 
£ using Eq. |3| A "direct" approach is to extract £ directly from the matrix elements using 



m Bd M 



M r . 



(6) 



With the notable exceptions 0, 0] previous studies of B° — B° mixingin lattice QCD 
have been done in the quenched approximation. It has been estimated 0,0 that even the 
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values of the ratios, which are less sensitive to systematic errors, are likely to have 10-30% 
quenching errors. In order to mitigate the effects of quenching errors, this calculation was 
done on lattices with two dynamical quarks. 

This paper is organized as follows. In Sec. El we discuss the lattice action used in our 
simulation. In Sec. II 111 we focus on the symmetries of the action. In Sec. II VI we review the 
connection between continuum QCD and continuum heavy-quark effective theory (HQET). 
In Sec. we discuss the matching between continuum HQET operators and their lattice 
counterparts. In Sec. I VII we review our methods of computing matrix elements on the lattice. 
We present the results of our calculation in Sec. IVHI 



II. ACTION 



It is convenient to separate the action used in our simulation into three different parts. 

S Sh eaV y ~\- Slight Sg aU g £l (7) 

where Sh eavy and Su g ht describe the heavy- and light-quark sectors and S gauge is the gauge 
part of the action. In this section we discuss the choices made for each part. Unless it is 
given explicitly we assume the lattice spacing a to be equal to 1 for the rest of this paper. 



A. Heavy Quark Action 



Heavy Quark Effective TheoryfHQET) is an expansion of the QCD Lagrangian in inverse 
powers of the heavy-quark mass 0, H,l|. This expansion is justified as long as the momenta 
involved are small compared to the heavy-quark mass. For B mesons the typical momentum 
transfer is of the order of Aqcd ~ 200MeV. Since this is in fact much smaller than the mass 
of the b-quark, we expect that the leading order of HQET provides a sufficiently accurate 
description of B physics phenomena. Thus, we will take the m& — > oo limit of HQET, also 
known as static-quark approximation. On the lattice the static-quark action is given by 

Sheavy = ^2Hx)V h(x), (8) 

x 



V^h(x) = h(x) - U>(x - 0)h(x - 0), 



(9) 



where Uq(x — 0) is the gauge link in temporal direction between sites x and x — 0. The 
static-quark field h satisfies 70/i = h. We will discuss the implications and the symmetries 
associated with this formulation in the next section. Here we just note that the static 
approximation describes exactly what its name suggests: a heavy quark fixed in space and 
propagating in the time direction only. The propagator is given explicitly by 



G Q (x,y) = 6(t x - t v )5(x - y) \ut(y,t y )Ut(y,t y + 1) . . . Ul(y,t x - 1) 



(10) 



The lattice static-quark action (jHJ) is 0(a)-improved as was pointed out in Ref. [l(|. Previous 
studies [ll, lH IS 13; which relied on the action in Eq. (JHJ), were plagued by large statistical 
errors inherent to this formulation. However, following the ideas developed in Ref. [15], we 
use a modified lattice formulation of the static-quark action, which circumvents the poor 
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signal-to-noise ratio while maintaining the key advantages. We replace the simple gauge link 
Uq(x) in the backward derivative in Eq. Q with an average of the 6 staples around that 
link 



V °( x ) = r 2 \u 3 (x)Uo(x+j)U](x + 6) + Uj(x - j)U {x - j)U 3 (x + 6 - j) 



Now our propagator for the heavy quark is given explicitly by 

G h (x, y) = 6(t x - t v )8(x - y) \vj (y, t y )V j (y, t y + l)... V j (y, t x - 1) 



;i2) 



B. Light Quark Action 

We will use the domain-wall fermion (DWF) formalism to describe light quarks on the 
lattice. The DWF formulation achieves improved chiral properties at finite lattice spacing 
via introduction of an extra dimension [la. Il7| . Though the DWF method does not provide 
exact chiral symmetry, the degree of chiral symmetry breaking and the resulting errors are 
under control and can be precisely quantified [liL 19j |. For a light quark with mass m q the 
DWF action is 



L s -1 



light 



x 



s=0 



v; 



-v; v+ + m 5 + p L dt - P R d^ | t.(x) 

+ S ^m q q(x)q(x). (13) 



Note that the fermion field ip s has an index s, which represents the fifth dimension. The 
fifth dimension extends from to L s — 1 and exact chiral symmetry is achieved in the 
L s — > oo limit. Gauge fields are confined to the four-dimensional boundaries at s = and 
L s — 1 and are present in the forward and backward covariant derivatives V^: 



Vti>i x ) = U^(x)i>(x + jX) - ip{x) 



(14) 



V„ ij)(x) = if>(x) - UUx - p)i/)(x - p). 



(15) 



Correspondingly df are derivatives in the fifth dimension with = 1. Pi = (1 — 75) /2 and 
Pr = (1 + 75)/2 are the projectors for the left- and right-handed spinors. The domain-wall 
height M 5 is a parameter of the theory, which we set equal to 1.8. Finally, q(x) represents 
the physical four-dimensional quark field constructed from the light modes at s = and 
L a -1 

q(x) = PlMx) + PriPl s -i(x), (16) 



q(x) = ip {x)P R + ipL 3 -i{x)P L . 



(17) 
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C. Gauge Action 



The simplest satisfactory formulation of SU(3) gauge fields on the lattice is given by the 
standard Wilson plaquette action 

6 P 

where (3 = 6/pq with g denoting the bare lattice coupling. Up is the path-ordered product 
of links around the lxl plaquette P. Recently, it has been realized that it is possible to 
construct lattice gauge actions that are superior to the naive plaquette action in terms of 
approxima ting the continuum gauge fields. We will use the doubly-blocked Wilson (DBW2) 
action 2(1 21L w hich was shown to have excellent chiral properties when combined with the 
DWF action 

Sgauge = ~f U 1 ~ 8c l) R ^l U P] + C l ReTr [^ > ( 19 ) 

where Ur is the path-ordered product of links around the 1x2 rectangle R. The parameter 
Ci was set to —1.4069 and (3 = 0.8 which corresponds to a -1 « 1.7 GeV. 



III. SYMMETRIES OF THE ACTION 



The advantages of the choices we made for our action in Eq. (J7J) become clear when we 
consider its large symmetry group. We focus on two symmetries particularly important for 
our calculation. 

A. Chiral symmetry 

SU \Nf) L ® SXJ(Nf)n chiral rotations are given by 

Ms) ^ U L {x)iP s {x), 0<s< L s /2 - 1, (20) 



1>s(x) -> U R {x)iP s {x), L s /2 <s<L s -l, (21) 

where Ul(x) and Ur(x) acting on fermion fields on the left- and right-hand halves of the 
five-dimensional lattice are given by 

U L (x) = e^ T \ (22) 



U R (x)=e-^ T \ (23) 

with T a being the generators of the SXJ(Nf) flavor group. For finite L s the symmetry due 
to transformations (}2T)j) and ()21j) is not exact. The slight breaking of the chiral symmetry 
manifests itself through the appearance of an additional mass term, in the effective 

Lagrangian. This residual mass term falls off rapidly with increasing L s [19]. For our choice 
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of L s = 12 we observe m res /m s = 0.03 or m res m 2.5 MeV. We will take into account the 
explicit breaking of chiral symmetry due to m res simply by shifting the values of our light- 
quark masses from m/ to rrif + m res . Note that the additional chiral symmetry breaking due 
to the appearance of the dimension- five operators is of the order am res AQ C D ~ 0.5 MeV. 
Since it is much smaller than other systematic errors, we will consider the effects of the 
dimension-five operators to be negligible. 



B. Heavy-quark spin symmetry 

The SU(2) heavy- quark spin symmetry (HQS) transformation is defined as 

V> fc (x) -> Vil> h {x), (24) 

with 

y = g-i&eyjto^ ^5) 

where fa is a parameter and = §[7i,7j]. The heavy-quark spin symmetry simply means 
that the heavy-light quark interactions are invariant under arbitrary changes of the spin of 
the heavy quark. Note also that due to the heavy-quark field equation 70/1 = h, the tensor 
density T^ v is not an independent bilinear. In particular, T - 7 = V J and T iJ = A . 



IV. MATRIX ELEMENTS IN THE STATIC LIMIT 

The expression for the hadronic matrix element in Eq. Q is a continuum QCD result. 
It is not obvious how to compute it directly on the lattice, since the inverse lattice spacing 
a -1 < mb for the current generation of computers. In order to make it computable we will 
make use of HQET. The relevant continuum operator, here normalized according to the 
MS(NDR) scheme, is 

0f^l A) = [W - 75)9] - 75)?] • (26) 

This operator is evaluated between parity-even states in Eq. (J2J). Therefore, only the parity- 
even part gives a nonzero contribution to the matrix element. The parity-even part 

O^KP = (67*9) (M + pfw) (h^q) , (27) 
can be written using the operator product expansion as following 

0^aa\^) = Z l (» b ^)0^f A M + Z 2 (^ b ^)O^M + 0{l/tit). (28) 

The leading dependence of Oyy^^ on the b-quark mass is absorbed into analytic coeffi- 
cients /x) and Z2 (//&,//). Effectively, encode the "short-distance" physics between 
scales //& « 5 GeV and fi ~ 1.7 GeV whereas Of QET (/z) are operators in HQET, which 
describe the "long-distance" physics below fi. 

O vv f AA = 2 g) (ti-\q) + 2 (/^y 7 5<7) (h^^q) , (29) 
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n HQET 
W SS+PP 



where 



h^\x) = e ±imv - x l±lflV \ (x). 



(30) 
(31) 



The full-theory operator in Eq. (|2*7|) can create two heavy quarks or annihilate two heavy 
antiquarks. Since the field h^ + ' cannot annihilate a heavy antiquark and h^~' cannot create 
a heavy quark we need a compensating factor of 2 in Eq. ()29|) and Eq. (|3U|). Note also that 
the operators (3 i H Q ET having no dependence on the mass of the b-quark are well suited for a 
static-limit lattice calculation. 

We quote the results for Z 12 obtained via perturbative calculations 0, 2^, 27] : 



Zi{m h ,n) 



a s (m b )~ 






25 ^ 



! _ u ^s}jnbl a a (jjt) -a s (m b ) j N 




a s (/x) 



Z 2 (m b ,n) = -8 



47T 



. a s (/i) 



" 25 



(32) 



(33) 



where we assumed the number of flavors rif = 4 for the evaluations of the anomalous di- 



mensions and beta-function coefficients. For J 



7<> 



we obtain 3.881. According 



to the latest PDG review Q A ( q ] cd = 217 MeV, which implies A^ D = 276 MeV for our 
calculation. Then a s {m b ) = 0.213 for m b = 4.3 GeV and a s (/i) = 0.309 for /i = 1.7 GeV. 
Thus, the numerical values for Z\ 2 are 



Z x = 0.880, Z 2 = -0.144. 



(34) 



Similarly, we can express the pseudoscalar meson decay constant fs in terms of the HQET 
matrix element 

/bv^ = Z (m b ,fM){0\hj -f 5 q\B), (35) 

where me is the meson mass and Z Q is known to the precision of the two-loop anomalous 
dimension calculation With the number of flavors equal to 4 the result is 



~a s (m b y 






25 ^ 



8 a s (m b ) a s {ji) - a s (m b ) ^ 
3 4tt 47r 



(36) 



with J 



7 (0) /3i 



^ = 0.948. Assuming other parameters to be the same as for Z 12 we get 



1.05. 



(37) 
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V. MATCHING CONTINUUM AND LATTICE OPERATORS 



In order to relate our lattice results to the physical values of observables we need to know 
how to match lattice operators to the effective theory operators defined in Sec. |TV| Such 
matching depends on the details of the action used in simulation. The matching coefficients 
calculated at one-loop level in perturbation theory or nonperturbatively are known for a 
limited number of actions. Unfortunately, a matching calculation for the combination of 
improved static quarks and DWF quarks has not yet been carried out. Therefore, in this 
paper we will do lattice-to-continuum matching at tree level. 

Previous studies involving static heavy quarks and Wilson fermions had significant correc- 
tions to the tree- level values of fs and Bb- Both perturbative 11, jy| and nonperturbative 
*32l 33l calculations indicate that for our a" 1 « 1.7 GeV there is 17 — 23% downward correc- 



tion to the tree-level value of fs- For the matrix element and Bb values there are downward 
20 — 37% corrections, depending on the way the perturbative calculation is performed, ac- 
cording to Ref. [3J, |36|. Since domain- wall fermions are akin to Wilson fermions, we 
expect that, apart from a peculiar wave function renormalization factor discussed below, the 
corrections to our tree- level values of fs and Bb will be similar. However, as will be argued 
below, the uncertainty in matching coefficients has no influence on the precision of /b s / fB d 
and has negligible influence on the values of BbJ Bs d and £. 



A. Jb matching 

For the case of the static heavy-quark and Wilson light-quark action the tree-level match- 
ing for the decay constant would be trivial. The matching coefficient Z q in Eq. would 
be 1. In the case of domain- wall quarks, the physical modes extend slightly into the fifth 
dimension; therefore, the overlap between the interpolating fields on the boundary and the 
physical modes is different from 1 even at tree level |37(. Fortunately, the quark field renor- 
malization factor, which includes this overlap effect, was calculated nonperturbatively in Ref 
j38|. With q representing the renormalized light-quark field we have 

q = Z^q 1 ^. (38) 

For our value of M 5 = 1.8 the Z q determined from the conserved axial current is 

Z]l 2 = 0.897(11). (39) 

Combining this result with the QCD-HQET matching described in Sec. II VI we obtain the 
final formula for fs 

f Bv ^ = Z Zl/ 2 ^a- 3 / 2 , (40) 

where $^ tt is the lattice counterpart of the continuum quantity /ba/wb- The numerical 
value of can be obtained via Eq. (jH2J). Note that by including Z q in Eq. (j4*Uj) we deviate 
from a strictly tree-level calculation. However, we stress that the Z q factor is specific to 
the domain-wall fermion formulation and contains a large tree-level piece. By including it 
we compensate for the fifth dimension and make our results more easily comparable to the 
studies done with Wilson fermions. Except for Z q there are no other renormalization factors, 
since all four-dimensional amplitudes are evaluated at tree level. 
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B. Four-fermion operator matching 



The tree-level matching of four-fermion operators is straightforward. There is no operator 
mixing and we only have to take into account the light-quark overlap factor. Thus 

= z?0 latt_ (41) 

Combining this with the QCD-HQET matching discussed in Sec. II VI we obtain 

n MS(NDR) _ 7 7 Watt \ 7 7 D latt (A9\ 

U VV+AA — Zj l Zj q U VV+AA + q U SS+PP' l 4Z J 

Note that Z q cancels in the Bb determination since 

/Rin M5(NDR) iR\ 

_ K^l^VV+AA \ D I _ ry 7 -2nl3tt i ry ry -2 olatt (An\ 

n B — |T2 2 ~~ ^ iZj ^VV+AA + ^2^0 n SS+PP' l 4J ) 

3 J B m B 

where 

olatt = (^|Ql att (0)l^) (u] 

°' §($^) 2 m B ' 1 ' 

Even though we match at tree level it is worthwhile to consider the effect of a one-loop 
matching. First, we observe that since Z 2 in Eq. (}3"3"j) is already order a s , the lattice matching 
at tree level is sufficient for 0^+p P . The operator Oyy^j^ needs to be matched at one- 
loop level. Because of chiral symmetry, Ovv+aa can mix only with Oss+pp and Ott- As 
discussed in Sec. lIIII Ott can be decomposed into Oss+pp and Ovv+aa- Therefore, Ovv+aa 
can mix with Oss+pp only. The perturbative calculations for Ovv+aa with Wilson quarks 
jjil 35[ show that it does not mix with Oss+pp- Thus, Ovv+aa should be renormalized 



by a single multiplicative constant. For a -1 rs 2 GeV the typical value of this constant is 
0.59. This factor would multiply Z\Z q in Eq. (|42|) and Z\ in Eq. (|43|) . With this correction 
our tree-level results for M. and Bb would decrease by 37%, but the results for the ratios 
M. s j 'Aid and BbJ BB d would change by less than 1%. Thus, for the purposes of the ratios 
determination the tree-level matching is sufficiently accurate. 

VI. COMPUTATION OF MATRIX ELEMENTS ON THE LATTICE 

Our method for computing the relevant matrix elements on the lattice is fairly standard. 
We calculate the following two- and three-point correlation functions 

C PC \t,t ) = ^(0|<(x,t)A c (t )t|0), (45) 

X 

C cc (t,t ) = (0\A^t)A^(t o y\0), (46) 

c 0i (h,t) = ]T<o|^(^o: at W)4Wlo>, (47) 

X 

with 

<(x,t) = /i(x,t)7o7 5 g(x,t), (48) 
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^o(*)= Yl h(x,t)jol5q(y,t), (49) 

(x,y)GC 

^o(*)= 9(x,t)7o75^(y,t), (50) 

(ar,j/)eC 

where C and P stand for cube-smeared and point interpolation operators used for B mesons. 
The smearing was done in a cube with side Lb = 9, since it showed the best isolation of the 
B-meson ground state. The sources were located at to = and t\ = 20. We used Coulomb 
gauge fixing for cube-smeared sources and sinks. 

The heavy-light two-point correlators for large t behave as 

C rC( t ,0)^A PC e- E \ (51) 

C cc( t ,0)^A cc e- Et . (52) 

We obtain A PC and A by simultaneous fit of C PC (t) and C cc (t) to an exponential function 
with E being the unphysical binding energy for the lightest pseudoscalar state. The quality 
of the ground state isolation can be monitored via plots of E(t) and Co, (20, t). As can be 
seen in FIG. [T]we have good ground state isolation starting at t ^ 6. Once A PC and A cc 
are determined we can calculate the matrix element of the lattice static-light axial current 

\/2A pc 

®T = m p *mB) = ^=. (53) 

The continuum HQET matrix element is then directly proportional to its lattice value 

(01^075915) = Z\' 2 ^\ (54) 

where Z q is a matching factor, which was discussed in Sec. fVl 

Similarly, the three-point function with sources for both heavy and light quarks located 
at and t\ behaves as 

C 0i (tut) -^-A cc (B\O^(0)\B)e~ Et \ (55) 

IVTiB 

The most elegant way to extract the matrix element would be to divide Coi(ti,t) by C cc (ti). 
Unfortunately, for t\ = 20 the smeared-sink two-point correlation function exhibits pro- 
hibitively large errors. We will circumvent this difficulty by using two two-point correlation 
functions at earlier times 

The raw data for Myv+AA an d Mss+pp are depicted in FIG. Eland FIG. El Finally, the 
lattice B parameter is calculated as follows 



p (+ + \- Cpj (*!»*) *i»<»0, rjlatt (Kt7 \ 

Roi {tut) = fc^MOc^M) B °< ■ (57) 



As can be seen in FIG. EJand FIG. EJthe errors for ifo. are smaller than for Mq^ since the 
point-sink two-point correlation functions are less noisy. 
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VII. NUMERICAL RESULTS 



We performed measurements on 16 3 x 32 lattices with two dynamical flavors of domain- 
wall quarks. The inverse lattice spacing determined from the p-meson mass corresponds 
to 1.691(53) GeV, which is consistent with the value of 1.688(21) GeV obtained from the 
Sommer scale. In order to make the chiral extrapolation all quantities were calculated at 
three values of the dynamical u and d quarks with the bare masses: 0.02, 0.03 and 0.04. The 
bare quark mass 0.0446(29) corresponds to the mass of the strange quark m s . The value for 
m res incorporating the slight breaking of the chiral symmetry was found to be 0.03m s . The 
sea and valence quark masses were set to be identical. The pseudoscalar pion mass for the 
three quark masses is 0.2910(24), 0.3568(25) and 0.4086(21) in lattice units. We performed 
calculations on 449, 425 and 488 configurations for rrif = 0.02, rrif — 0.03 and rrif = 0.04 
respectively. Configurations were chosen at random from 5,000 thermalized hybrid Monte 
Carlo trajectories separated by 10 trajectories. Statistical errors were computed using the 
jackknife method with block size equal to 5 configurations. Static-quark propagators were 
computed according to Eq. (|11|) and Eq. (|12j). For the light quarks, all Dirac matrix inver- 
sions were performed using the conjugate gradient algorithm with a stopping condition of 
10~ 8 . We used periodic boundary conditions in the spatial directions and antiperiodic in the 
time direction. More information on m res , conjugate gradient algorithm, a -1 determination 
and lattices generation can be found in Ref. [39j. 

Now we quote our results. The numerical values for the lattice $^ tt matrix element 
defined in Eq. (J53|) are given in Tab. 0for each value of the bare quark mass. The linear ex- 
trapolation to the chiral limit gives = 0.274(8) (22). Here and in all following results, the 
number in first parenthesis always represents statistical error, the numbers in other paren- 
theses represent systematic errors which will be discussed in detail at the end of this section. 
The value at the physical strange quark mass is 

= 0.356(8)(29)(5). Matching the lattice 
and continuum operators using Eq. flUJ) we obtain f Bi ^rn^ d = 0.568(17) (50) (^ 31 )GeV 3/2 , 
which corresponds to fs d = 247(7) (22) (^57) MeV where we used rriB d = 5279 MeV (see Ref. 
[Hf). We plotted /sy / ms in physical units as a function of mass in FIG. El 

We quoted the fs d value solely for the purpose of cross-checking our results with previous 
calculations in the static limit. Duncan et al. [11[ report fs d = 188(23) (^29) MeV with the 
lattice-to-continuum matching factor Z = 0.77 at a -1 = 1.78(9) GeV. We compensate for Z 
to make a comparison at tree level and obtain the central value of fs d = 244 MeV, which is 
quite similar to our result. 

Our main focus is the determination of the SU(3) breaking ratios which we can determine 
much more precisely. For the ratio of the decay constants we obtain 



where we assumed mB 3 /mB d = 1.017 (see Ref. |28(). Our value is somewhat larger than the 
conventional value fsj fB d = 1-15. Interestingly, McNeile and Michael [HI also using an 
improved lattice formulation of static approximation and unquenched clover action obtained 



The numerical values of Mo { defined in Eq. (I56j) are given in Tab. ITTlin lattice units. The 
physical hadronic matrix element Ai q can be obtained from Mo 4 via Eq. (|4*2*|) . We show the 
dependence of Aiqin^ 1 on the mass of the light quark in FIG. [7| The linear extrapolation 




(58) 



f B Jf Bd = 1.38(13)(8). 
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to the chiral limit gives Mdm^ 1 = 0.727(96) (65) (^ 69 ) GeV 3 and to the mass of the strange 
quark M. s m]f = 1.73(9)(15)(+g 4 )(7) GeV 3 . Thus, for the ratio of matrix elements we obtain 

£^ = 2.42(32) (21) (10). (59) 

To our knowledge, this is the first unquenched calculation that obtained the ratio of matrix 
elements directly. To put our result in context we quote results from a previous quenched 



study. Bernard, Blum and Soni [4l| . using extrapolated Wilson action found M. s /M.d = 
1.76(10) and 2.21(42) for constant and linear continuum extrapolations, respectively. Using 
Eq. © we calculate £ directly from the ratio of the matrix elements 

e= 1.53(10)(7)(3). (60) 

Finally, the values for Bq^ obtained using Eq. (J57|) can be found in Tab. 11111 The dependence 
of Bb obtained using Eq. (|43J) on the mass is depicted in FIG. |H1 Again using the linear 
extrapolation we obtain B Bd = 0.812(48) (67) (tjj 00 ) and B Ba = 0.864(28) (71) (^ 20 ) (3). Thus 

Bb. 



1.06(6)(3)(1). (61) 

ratio with the ratio 

rect" way via Eq. © 



B Bd 

Combining the above ratio with the ratio of the decay constants we calculate £ in an "indi- 

i 11 T71 TPlv 

e= 1.33(8)(6)(2). (62) 

As can be seen in FIG. El and FIG. El the separate linear extrapolations for Jb^^ib and 
Bb fit data points better than the linear fit for AAm^ in FIG. Therefore, the "indirect" 
method seems to be more reliable than the "direct" method for the case of linear chiral 
extrapolation. Thus, we adopt as our final answer for £ the value in Eq. (|62p. 

Now let us discuss errors. The first error is statistical and is easily quantified using the 
jackknife method. The number in second parenthesis is the finite volume error and the finite 
lattice spacing error added in quadrature. As argued in Ref. |42| finite volume effects for the 
above quantities are small (^ 2%) as long as the calculation is unquenched. Finite volume 
errors are being estimated using x?T which we believe should be accurate for mesons in our 
volume. The finite lattice spacing errors present a larger concern. Though the action we 
use is 0(a)-improved, the operators themselves are not. Therefore, 0(a) errors are possible 
for the individual quantities. A continuum extrapolation reported in Ref. [43| suggests that 
the error can be as large as 8%. Note that these estimates of finite lattice spacing errors 
themselves contain quenching errors which could be large. For the individual quantities the 
biggest uncertainty by far comes from the corrections to the tree-level matching between 
the continuum and lattice operators. As discussed in Sec. [V] we expect these corrections to 
be similar to the ones in the case of static heavy quarks and Wilson fermions. We adopt 
the largest published corrections [ill . 34| to tree-level as our error estimates. Thus, the 



uncertainty associated with the lattice-to-continuum matching is 23% for f B and 37% for 
the matrix elements and Bb- This matching error is quoted in the third parenthesis for 
fs-i -M and Bb- For the quantities involving the strange quark, there is an uncertainty 
associated with the value of the strange quark mass. It is quoted in the third parenthesis 
for the ratios. It is quoted in the fourth parenthesis for /^ s , M-s and Bb s - 

Finally, let us discuss the uncertainty associated with chiral extrapolation. The chiral 
correction to Jb predicted via is well known j^, 0|. Here, for the case of Nf = 2 we 
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have 



Ib — fo 



3(1 + 3^) m PS m PS 
4 (4tt/) 2 /x 2 



(63) 



where g is the heavy-meson coupling to pions. The same expression with —(1 — 3g 2 )/2 
instead of 3(1 + 3g 2 )/4 is predicted for Bb- A possible approach to chiral extrapolation 
is to fit the lattice data to Eq. (J63|) and then to extrapolate to small values of pion mass. 
However, as was pointed out in Ref. |45| such an approach is misleading since Eq. (JB13|) is not 
valid at large values of pion mass (nips > 500 MeV). On the other hand, we realize that the 
linear chiral extrapolation used in our calculation is not satisfactory either, since it ignores 



legitimate corrections coming from long-distance physics [46, 47]. Clearly, the issue can be 
resolved only when more data points at lower values of quark masses become available. In 
our case, an attempt to quantify the chiral extrapolation error via xPT would border on 
speculation, since our data is consistent with the linear extrapolation and it is outside the 
region of validity of xPT. Therefore, we do not include the uncertainty associated with the 
chiral extrapolation in our estimation of systematic error. 

Our calculation was performed in the limit of the infinite heavy-quark mass. In principle, 
the finiteness of the heavy-quark mass can be taken into account by including 1/mg correc- 
tions. The naive expectation is that the corrections to the individual quantities should be 
of the order Aqcd/ttib- Then the corrections to the values of the ratios should not exceed 
Aqcd/itt-b times the difference between the value of the ratio and 1. For example, in the 
case of £, the corrections can be of the order (£ — 1)Aqcd/™b. Thus, if l/ms corrections are 
taken into account, the change in the value of Bp ratio should not be noticeable, whereas 
the values of the Jb ratio and £ could be modified by 2% — 4%. 

In summary, our final answers for SU(3) flavor breaking ratios are: fB s /f'B d — l-29(4)(6), 
Bbs/ Bb a = 1.06(6)(4), £ = 1.33(8)(8) where the number in the first parenthesis is statistical 
error and the number in the second parenthesis is the arithmetic sum of systematic errors. 
Note that these systematic errors do not include the uncertainty associated with chiral 
extrapolation and possible I/tub corrections. 



VIII. CONCLUSIONS 

In this paper we have done an unquenched two-flavor calculation of SU(3) flavor breaking 
ratios in B° — B° mixing. A chirally improved, domain-wall fermion action was used for 
the light quarks. An improved lattice formulation of the static action was used for the 
heavy quarks. Two flavors of dynamical quarks were simulated at three values of quark 
mass with the lowest value equal to half of the physical strange quark mass. Lattice-to- 
continuum matching of operators was done at tree level, which resulted in large uncertainties 
for individual quantities, but had negligible influence on the value of the ratios due to the 
large symmetry groups of the action. Our results suggest that the conventional values for 
phenomenologically relevant parameters fB s /fB d and £ may be too low. 
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TABLE I: The values of defined in Eq. (|53|) for different light quark masses. Statistical errors 
are given in parentheses. 



m/ + m res 


fftlatt 


0.0214 
0.0314 
0.0414 


0.314(4) 
0.330(4) 
0.351(4) 


TABLE II: The values of lattice matrix elements defined in Eq. I)56|) for different light quark 
masses. Statistical errors are given in parentheses. 


mf + m res M V v+AA 


Mss+pp 


0.0214 0.326(13) 
0.0314 0.366(11) 
0.0414 0.446(12) 


-0.206(7) 
-0.235(7) 
-0.276(7) 


TABLE III: The values of B l gf defined in Eq. (JSJ) and obtained 
quark masses. Statistical errors are given in parentheses. 


via Eq. (|57|) for different light 


m f m res By y +J 4^4 


73 latt 
n SS+PP 


0.0214 0.960(28) 
0.0314 0.950(18) 
0.0414 0.987(16) 


-0.605(15) 
-0.613(10) 
-0.611(9) 
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FIG. 1: The binding energy for the lightest pseudoscalar state E(t) = ^ g £S( f+1 ) an d the three- 
point correlation function Cyv+AA (20, t) as functions of time. Both plots are for the bare light 
quark mass rrif = 0.04. Error bars represent statistical errors. 
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FIG. 2: A plot showing the raw data for M Ovv+AA (20,t) defined in Eq. (J5H and Eq. (j2HJ) for 
different values of light quark masses. Error bars represent statistical errors. Plateaus indicate the 
weighted averages for 6 ^ t ^ 14. Their values are quoted in Tab. ITTl 
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FIG. 3: A plot showing the raw data for M Oss+PP (20,t) defined in Eq. (JM) and Eq. (jHOJ for 
different values of light quark masses. Error bars represent statistical errors. Plateaus indicate the 
weighted averages for 6 ^ t ^ 14. Their values are quoted in Tab. ITT1 
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FIG. 4: A plot showing the raw data for R Ovv+AA (20, t) defined in Eq. (jnZJ) and Eq. (EU) for 
different values of light quark masses. Error bars represent statistical errors. Plateaus indicate the 
weighted averages for 6 ^ t ^ 14. Their values are quoted in Tab. Hill 
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FIG. 5: A plot showing the raw data for Ro ss+PP (20, t) defined in Eq. (|STj) and Eq. (j50J) for 
different values of light quark masses. Error bars represent statistical errors. Plateaus indicate the 
weighted averages for 6 ^ t ^ 14. Their values are quoted in Tab. Illll 
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FIG. 7: Chiral extrapolation of A4 q m~^ obtained using Eq. I|42|). Error bars represent statistical 
errors. 
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